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We propose to search for a tetraquark candidate X{udss) in pn ~* AAX (udss) AAK'^ 
or AAKK* . The existence of tetraquark state X{udss) with = 0^, 1~ or 1"*" was predicted 
in the hterature based on specific diquark efi'ective degrees of freedom inside hadrons. In order 
to understand the underlying dynamics for exotic hadrons, a search for the tetraquark X{udss) 
is strongly recommended. The proposed reaction involves two A production, of which the narrow 
widths make it a great advantage in the analysis of the final state missing mass spectrum. We make 
an estimate of the production rate of X{udss) in an effective Lagrangian theory and find that for 
= 1~ the sample events of ~ 2200 nb~^ will be able to identify X{udss) with five standard 
deviations at a width of 10 MeV to K+K° near threshold. For = 1+ with a width of 20 MeV to 
KK* , the sample events of ~ 130n6~'^ will be needed. Large production cross sections are expected 
in a kinematic region beyond the threshold. We emphasize the advantage of low background in this 
transition channel, and in the meantime caution the large uncertainties in the present estimate due 
to lack of knowledge about the X{udss) state. Implications for its heavy-fiavored partners qqcc and 
qqbb are briefiy discussed. 



PACS numbers: 13.25.-k, 13.75.Ev 



I. INTRODUCTION 



One puzzling question in low-energy QCD is the apparent absence of "exotic states" in hadron spectrum. 
Here, "exotic states" are referred to baryons beyond qqq or mesons beyond qq which are allowed by 
QCD symmetry but different from conventional constituent quark model classifications. So far, although 
there are some candidates for such exotic states, none of those has been indisputably established in 
experiment 

For those exotic states, if they have quantum numbers that cannot be constructed by conventional 
qqq or qq, determination of their exotic quantum numbers is a direct way of confirming their existence. 
However, for those with conventional quantum numbers, evidence for their existence can be contaminated 
by conventional quark model configurations due to possible configuration mixings. Their determination 
hence is a challenge for both experiment and theory. In 2003, experimental evidence for a pentaquark 
0"*" {uudds) was reported by LEPS Collaboration at SPring-8 2], which immediately initiated tremen- 
dous interests and activities in both experiment and theory. With a flavor configuration of uudds, i.e. 
strangeness S = +1, confirmation of such a state would be a direct evidence for the existence of exotic 
states. 

Although the final confirmation of this pentaquark state still needs further experimental studies, it 
sparkled many novel ideals in the understanding of quark-quark interactions and effective degrees of 
freedom within hadrons, among which the diquark effective degrees of freedom have attracted a lot of 
attention [3, 0, We shall not review those progresses here since there have been a lot of detailed 
discussions in the literature. We only point out that as a follow-up of the pentaquark scenario it will give 
rise to specific properties of effective diquark degrees of freedom, hence will predict existence of other 
multiquark exotics such as tetraquark states due to the attractive forces between the diquark clusters fo'] . 
It is interesting to realize that such diquark effective degrees of freedom may lead to a formation of 
broad tetraquarks even though a narrow pentaquark presumably does not exist. Studying their possible 
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manifestations in experiment forms our motivation of this work. In particular, we shall identify a unique 
channel for a tetraquark candidate production, and investigate possibilities of establishing or eliminating 
it in future experimental measurements. 

A promising tetraquark candidate with the valence quark content udss has been broadly studied 
recently 0, 0,11, i, [iO, [ll, while some relevant earlier work has been done by Jaffe decades ago [;3|. 
Burns, Close and Dudek [3| suggested that if the diquark degrees of freedom such as the ones proposed 
for the pentaquark (uudds) [3, IE] exist, there should also exist an isoscalar tetraquark udss with 
= 1", mass about 1.6 GeV, and decaying into K^K+ with a width around 0(10 100) MeV. They 
used a novel configuration to turn a "bad diquark" into a "good diquark" . Karliner and Lipkin Q thought 
there would be an isoscalar udss meson with = 0+. As it cannot couple to K^K~^ or KKn considering 
the generalized Bose statistics and parity conservation, the lowest decay mode would be four-body KKttt: 
channel which may result in a narrow width and make the state detectable. In contrast, Kanada-En'yo 
et al. argued that the udss meson with = 1+ may be a stable and low-lying state, and decay into 
KK* via 5- wave |8|]. 

Motivated by the above studies, we propose to search for this tetraquark candidate X{udss) in pn 
AAX (udss) which is illustrated in Fig. [T] Although there are still large uncertainties with the model 
predictions, the cross section may be small and the width of X{udss) may be broad, the background 
which contains two narrow A's, is so clean that it may be possible to extract resonance signals in the 
missing mass spectrum of two A final states. 

II. THE MODEL 

In this section we shall study tetraquarks of = 1^ and 1+ as two most possible candidates with 
effective Lagrangians. In principle, the tetraquark state can be produced via several typical transitions 
as illustrated in Figs. [T] and [2] which can be labeled by the Mandalstam variables, i.e. t, u and s-channel. 
In Fig. [T] (t-channel), two meson exchanges are required. In Fig. [U^a) (w-channel) a meson and a S state 
will mediate the transition while in Fig. EJb) (s-channel) a meson and a pentaquark state 9 exchange are 
needed. Since we have no applicable information on the u and s-channel transitions, we do not discuss 
their contributions in this work. We shall focus on the i-channel transition and the details of the model 
are given as follows. 

A. Production of X(udss) with = 1~ 

For the process illustrated in Fig. [U the following effective Lagrangians are adopted for the production 
of A (udss) with = 1": 

Cknk = -igKNANj^AK + H.c. (1) 

-C/f-WA = -9K'NaN (^^^AK*'' - !l^-^a,,,Ad''K*^^^ + H.c. (2) 

Ckkx = igKKxX^" (d^K^K- - K^d^^K-^ + H.c. (3) 
'^''^dc^Xpd^KIK + H.c. (4) 

Then we adopt the following values for the KNA and K*NA couplings from Refs. [l2,[i3| which are also 
adopted by Ref. 14]: 

gK-NA = -4.26, KK'NA = 2.66, 

9kna = (1 + '2F)g^NN = -13.24, (5) 
where g^jvA obtained by using SU(3) flavor symmetry relation. 
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The coupling constant qkkx is determined by the decay width of X which is predicted by Ref. Q: 
V{X K^K^) ^ 0(10 ~ 100) MeV. With the effective Lagrangian mentioned previously, we have: 



^ ^ 487r ' ^ ' 



9kkx {m\ - '^K?'^ 
A 

which leads to 

gKKX = 1-4, when Y[X -> K^K^) ^ 10 MeV. (7) 
We also assume that Y{X KK*) ^ 0(10 ~ 100) MeV, with the similar process, we have: 

9kk*x [{^\ - i'^K' + mKf){m]^ ~ [mK- - mK)2)]i/2 



T{X KK*) = 



487r 2m\ 



X [{m,j^ + ml, - mlf - Am\ml,], (8) 

and 

giiii'X = 2.6, when T{X -> KK*) = 10 MeV. (9) 

For the vertices KNA and K*NA, we introduce the covariant monopole form factor at each interaction 
vertex as the intermediate meson may be off-shell. 



A2 - M, 



where A — Akna, Ak'NA here, and M^x are the corresponding masses of the exchanged mesons, uik or 
rriK'- The following commonly used values are adopted jTlj: 

Akna = 1.1 GeV, A^-iVA - 1-0 GeV. (11) 

So far there are no available data for the determination of the KKX and KK*X form factors. We 
hence introduce a dipole form factor for the meson coupling vertices, i.e. 

ft(.')^(^^)'. (12) 

where Ax is the cut-off energies for KKX and KK*X, and M^x denotes the corresponding masses of 
the exchanged mesons, rriK or tuk' , respectively. The Ax-dependence of the total cross section will be 
illustrated later. 

In Fig. [1] some of the kinematic parameters are defined as follows: 

qi ^ ri- si, q2^r2- S2 

Wi = ^1-52, W2^Si-r2. (13) 



The corresponding transition matrix elements are then given as follows 

(gi - g2) ■ ex 
' {<ll - ml){ql - ml) 

(wi ~ W2) ■ ex 
' (wf - ml){wj - ml) 



Tifi = gKKxgl:NAu{si)l5u{ri)u{s2)l5u{r2)— 2w 2 _ — 2\PK{ql)FK{ql)FKKx{ql) 

gKKxgj(NA'^is2)-f5u{ri)u{si)-f5u{r2) ^ Fk (wf) Fk (wj) Fk K X (wj) , (14) 



T2^f = Tkna (si,ri)r^.^A (S2,r2)r^^^.^ (^0,92) / 2 ^^2S^2^'^^ 2 s FKiqDFK'iqDFKK'xiql) 



-Fkjva (s2, ri) r^. (si, r2) ^kk-x i^o, W2) , ^ ^^^2\("^2 2 ^ Fk{wI)Fk' {wI)Fkk'x{wI), 

(15) 
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where 



r_R-ArA isi,ri) = -igxNAU (si) 75M (ri) (16) 

r^.ArA(.2,r2) = gK.^A^(.2) (7^ - ^^^^^^f^^^^) " (^2) (17) 
r^ifJC (*0, 92) — gKK'X£' a. 927 (18) 

PM^g.s-'^-^ ■ (19) 

The transition matrix element T3 for Fig. [T](c) can be obtained by making a momentum substitution in 
T2/i, i.e. ri ^ r2 and si ^ S2- 



B. Production of X{udss) with = 1+ 

As studied in ,8,]j axial vector meson udss {J^ ~ 1^) is a good tetraquark candidate. Within a flux- 
tube model, this state might appear around 1.4 GeV with a width of 0(20 ~ 80) MeV. We apply these 
quantities as an input to give an estimation of the cross section. The process is similar to the previous 
section except for the effective Lagrangian for the coupling of KK*X: 

^KK'X = gKK'X {d°'K*f^do,KX0 - d°'K*'^dpKX^) + H.c. (20) 

The meson X{udss) with = 1+ can be produced via ifi^T* -exchange but not ifif-exchange in this 
process. Thus we only need to consider Fig.[T](b) and (c). The following is the transition amplitude: 

T2./ = TkNA (Sl, ri) r^.^TA (52,^2) r^*A'*X (^0, 92) / 2 ^^''2U^2''^^ 2 ^ FKiqDPK* {ql)FKK'x{ql) 
-TkNA (s2,7-l)r^.ArA isi,r2)TKK'X (^0,^2) , ^ ^^^2\("^2 ^''^ 2 ^ ^^(^1 )fg» (wi)-FKK*x(w|), 

(21) 

where 

^KK'X = 9KK'xiqi ■ ^2/' " q^l), (22) 

and the other terms are the same as the previous section. The coupling constant gxK'X is determined 
by the width: 

T{X ^ KK*) = [2{pK-PK'? + ml,{ml + Ipkc^H] , (23) 

which leads to 

gKK-x = 7.5, when T{X -> KK*) = 20 MeV. (24) 

It should be cautioned that the adopted couplings still bare large uncertainties due to our limited 
knowledge on the tetraquark states. This will consequently bring uncertainties to the estimated produc- 
tion cross sections. However, we emphasize that our strategy here is to single out the pn KA.K^K^ 
and AAKK* channel which are advantageous for detecting the tetraquark states in experiment. Our 
calculation is to provide a reasonable estimate of the feasibility for future experiments. 
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C. Numerical Results 

Experimental data for pn — > AAK'^K'^ and AAKK* are not available so far. Interestingly, there are 
a few experiments in 1980's on the inclusive reaction pp — > 2A + anything [l^ and exclusive reaction 
pp ^ 2A + 2K~^ [H, UH motivated by the search for dibaryon state formed by two A. The total cross 
section of the inclusive reaction is about 19 ± 10 with the beam energy at £'„ — 32.1 GeV {W = 7.87 
GeV), and that of the exclusive reaction has an upper hmit about 460 nb at £'„ = 7.8 GeV {W = 4.05 
GeV). These data still possess large uncertainties and need to be improved. But at this moment, they 
can serve as a guidance for constraining our parameter space, and give an estimate of the cross sections 
for X tetraquark productions. 

Although the values for coupling constants can be determined indirectly by other experimental data 
such as strangeness productions and theoretical model prediction for X K^K'^, we still lack infor- 
mation about the choice of the cut-off energies in the form factors. Therefore, a thorough investigation 
of the parameter space is necessary. We first fix the parameters Akna, Ak'na as in Eq. pip and set 
Akkx = Akk'X = Ax = 1.2 GeV, which is a commonly adopted value for meson-meson interaction 
form factors. This allows us to obtain the total cross sections for different X states. We then examine the 
model sensitivities to Ax by extracting the Ax dependence of the total cross sections at a given energy. 

In Fig. [21 the energy dependence of the total cross section for pn — > AAX is presented for X with 
= 1^. It shows that the contribution is dominated by i^if-exchange, while KK*-eyic\\a.Tige is small. 
The total cross section also exhibits enhancement above threshold at ~ 7.0 GeV, and then dies out 
with the increasing W . The peak value is about 0.115 which is much smaller than the cross section 
for pp AA + anything [ITI. T his value is also below the upper limit of the exclusive pp ^ 2A -\- 2K^ 
cross section near threshold |18| . As shown by the solid curve in Fig. [31 the cross section near threshold 
at En — 7.8 GeV {W — 4.05 GeV) is about 2.3 nb. This value can be seen more clearly in Fig. [HJa) at 
Ax = 1.2 GeV. 

The sensitivity of the total cross section to the cut-off energy Ax is examined in Fig. [H The range 
of Ax — 1.0 ^ 1.5 GeV is a commonly accepted one in the literature. As shown by the solid curve, 
the cross section varies between 0.05 and 0.3 /x5 in terms of Ax, which suggests some sensitivities to the 
cut-off energies. The Ax-dependences of the exclusive KK and KK* cross sections are also presented 
by the dashed and dotted curves, respectively, and a relatively sensitive behavior of the KK* exchange 
is found. However, since contributions from this transition are rather small, the overall behavior of the 
cross sections is dominated by the KK exchange. Although the sensitivity to the cut-off energy brings 
some uncertainties to the estimate, it does not prevent us from drawing some preliminary conclusions on 
the production rate for X. 

Figure [5l shows the results for X production with = 1+. The underlying transition is via the KK* 
exchanges and the cross section increases with the energies. A fast rise appears near threshold, and then 
the cross section becomes rather flat. This behavior is different from the case of = 1~, where an 
obvious enhancement appears above threshold. At < 10 GeV, the total cross section is less than 0.21 
//6. Then, it slowly increases to 0.25 fib a.t W = 20 GeV. Nevertheless, the cross sections for 1+ are 
relatively larger than those for 1~ over a wide range of W. The different behavior of the near threshold 
cross sections between 1~ and 1"*" makes it an interesting place for looking for the X state in experiment. 

The Ax-dependence of the total cross section is also investigated and the results are displayed in Fig.[H] 
The cross section turns to be more sensitive to Ax than the case of = 1^. This will increase the 
uncertainties of the model predictions. However, in terms of gaining a rough idea about the production 
rate of X and providing a guidance of the future experimental plan, the range of the uncertainties can 
still be regarded as acceptable. 

Taking the experimental data for pp ^ 2A + 2K~^ [l3| as a guidance, we can estimate the number of 
sample events for establishing the X{udss). It is reasonable to assume that the cross section from the 
background contributions is atkg = 460 nb {W — 4.05GeV), while the signal cross section is ax — 2.3 nb 
[J^ = 1~) or ax — 9.5 nb [J^ = 1+). With a luminosity of L and event-collecting time t, we require 
the tetraquark signal of five standard deviations (5(t): 



(25) 
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where Nx and N^kg are the sample events for X{udss) and the background, respectively, and they are 
given by 

Nx = Lxtxax , (26) 

Nbkg EE LxtX (Jbkg ■ (27) 

This leads to: 

Nbkg > 1.0 X 10^ or Nx > 5.0 x 10^ or i x t > 2174 nb^\ for = 1" , (28) 

and Nbkg > 5.8 x lO'^ or Nx > 1.2 x 10^ or L x t > 127 for = 1+ . (29) 

Note that the above estimations are in the near-threshold region with W = 4.05 GeV instead of the 
kinematics with the largest cross sections, e.g. ~ 7 GeV. The reason is that we would like to compare 
our predictions with the only relevant experimental data from Ref. [l8l |. of which the inclusive cross 
section provide a rough check for the self-consistency of our model. It would be more interesting to 
look at the energy region oiW = 7 GeV where the cross sections are predicted to have a maximum for 
= 1~ and be sizeable for = 1"*". Given the experimental availability in the future, larger cross 
sections aXW = 7 GeV will make the detection of X much easier there. 

III. SUMMARY 

In this work we studied the production of possible tetraquark candidates X{udss) in pn — > A AX. 
Difficulty in the search for such an exotic state lies on the generally-large background in its productions, 
while as here we propose that pn AAX (udss) KKK^K^ and KKKK* are rather clean and ideal 
for looking for its signals. 

The numerical results for the cross sections were presented with a reasonable consideration of the 
parameter spaces. For X{udss) of = 1~, an obvious enhancement were observed above threshold, 
while for = 1+, no enhancement was seen. In both cases, typical values of hundreds of nb were found 
for the total cross section above the threshold region. 

We adopted the experimental upper hmit of 460 nb at E„ = 7.8 GeV {W = 4.05 GeV) [ll|, and 
estimate the signal-background rate. It shows that in order to observe the signal of X{udss) at the 
standard deviations of 5a, the sample events of the background must be larger than 1.0 x 10^ for = 1~ 
(or 5.8 X IC* for = 1+), or the product of luminosity and experiment time should be larger than about 
2174 nb~^ for — 1^ (or 127 nb^^ for = 1+). It should be noted that the estimates were made near 
threshold where the cross sections are not the maximum. This is due to the consideration of adopting the 
only available experimental information of Ref. [isj to estimate the signal-background rate. To search for 
the signal of X{udss), the ideal kinematic region should be around ~ 7 GeV, where an enhancement of 
the cross section was predicted for = l~ and the cross section also turned to be sizeable for — 1+. 
Meanwhile, we caution that although the theoretical study of the tetraquark properties is more essentially 
based on the diquark degrees of freedom, our knowledge on the pentaquark should also have influence 
on the estimate of the tetraquark decay widths. This will bring model-dependence to the theoretical 
predictions. Taking into account the form factors and couplings, the uncertainties of the results can be 
as large as two orders of magnitude. 

The similar experimental scheme could also be used to search for other exotic states such as X{qqcc) 
in the process NN — ^ AcAcX or X{qqbb) in the process NN AbAbX and so on. Since we lack experi- 
mental constraints on the DNAc {BNAb) and D*NAc {B*NAb) couplings, we cannot make quantitative 
estimations for those tetraquark productions. But we stress the advantages of such a reaction process for 
the search for tetraquark species of X{qqQQ). Experimental search for such exotics in hadron collider 
should be able to provide deeper insights into the properties of strong QCD dynamics. 
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FIG. 1: The tetraquark candidate X{udss) production via t-channel transitions in pn AAX. 
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FIG. 2: The tetraquark candidate X{udss) production via (a) u-channel and (b) s-channel transitions in pn 
A AX. 
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FIG. 3: Energy dependence of total cross sections for pn AAX with = 1 for the X. We have assumed 
r{X K+K°) = r{X KK") = 10 MeV Q. The cut-off energy is set as Akkx = Akk-x = Ax = 1.2 GeV. 
The dashed line is the contribution of ifA'-exchange, dotted line is that of A"* -exchange, and solid line is the 
total contribution. 
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FIG. 5: Energy dependence of total cross sections for pn — > AAX with = 1+ for the X. We have assumed 
r{X KK") = 20 MeV d- The cut-off energy is set as Ax = 1.2 GeV. 



